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Abstract. We present the results of analysis of constraints on cosmological parameters from cosmic microwave 
background (CMB) alone and in combination with galaxy cluster baryon fraction assuming inflation-generated 
adiabatic scalar fluctuations. The CMB constraints are obtained using our likelihood approximation method 
(Bartlett et al., 2000, Douspis et al., 2001). In the present analysis we use the new data coming from MAXIMA 
and BOOMERanG balloon borne experiments, the first results of the DASI interferometer together with the 
COBE/DMR data. The quality of these independent data sets implies that the Ct are rather well known, and 
allow reliable constraints. We found that the constraints in the SI — Ho plane are very tightened, favouring a flat 
Universe, that the index of the primordial fluctuations is very close to one, that the primordial baryon density is 
now in good agreement with primordial nucleosynthesis. Nevertheless degeneracies between several parameters still 
exist, and for instance the constraint on the cosmological constant or the Hubble constant are very weak, preferred 
values being low. A way to break these degeneracies is to "cross-constrain" the parameters by combining with 
constraints from other independent data. We use the baryon fraction determination from X-ray clusters of galaxies 
as an additional constraint and show that the combined analysis leads to strong constraints on all cosmological 
parameters. Using a high baryon fraction (~ 15% for h = 0.5) we found a rather low Hubble constant, values 
around 80 km/s/Mpc being ruled out. Using a recent and low baryon fraction estimation (~ 10% for h — 0.5) 
we found a preferred model with a low Hubble constant and a high density content (Q m ), an Einstein-de Sitter 
model being only weakly ruled out. 
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1. Introduction 

The determination of cosmological parameters is a cen- 
tral goal of modern cosmology. The measurements of the 
angular spectrum of the fluctuations in the CMB is one 
of the most promising techniques. Indeed the first mea- 
surements of fluctuations at the degree scales (Saskatoon, 
Netterfield et al., 1997) has allowed to put the first con- 
straints on cosmological parameters (Lineweaver et al. 
1997; Hancock et al 1997). Probably the most spectac- 
ular result was that the open models could be reasonably 
excluded, while flat models were favoured (Lineweaver & 
Barbosa; 1998b), essentially because of the location of 
the first peak. This result has been brilliantly confirmed 
by the first BOOMERanG and MAXIMA measurements 
(Hanany et al., 2000, de Bernardis et al., 2000). The data 
are now becoming of high quality and the detailed fea- 
tures expected in so-called inflationary models are now 
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becoming apparent (the famous peaks) from indepen- 
dent measurements (BOOMERanG, MAXIMA, DASI). 
Alternative scenarios like topological defects models are 
almost entirely ruled out. The possibility of a partial con- 
tribution remains open (Bouchet et al., 2000) and will be 
difficult to entirely rule out, although inflationary models 
are clearly preferred. In this paper we analyse the con- 
straints that can be set from the recent measurements 
obtained by BOOMERanG, MAXIMA, DASI combined 
with the COBE/DMR results. 

However, it has been realised that combination of dif- 
ferent methods to constraint cosmological parameters was 
very welcome, not only because of the gain in precision, 
but mainly because most of cosmological tests present de- 
generacies and may constraint rather well only specific 
combinations of cosmological parameters but not each pa- 
rameter individually. Combinations of constraints or mea- 
surements, often implemented as prior on some param- 
eters in likelihood analysis, are currently used to obtain 
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Fig. 1. The power plane: measured flat-band power 
estimates (COBE, BOOMERanG, MAXIMA and 
DASI) and our best fit model (by approximate 
maximum likelihood): (H , J7tot, A , Qbh 2 , n, Q) = 
(30 km/s/Mpc,1.3, 0.1, 0.019, 0.91, 23.0 /xK). 

accurate values for the cosmological parameters. In this 
paper, we use as a first additional constraint the observed 
baryon fraction in clusters. This quantity has the advan- 
tage of not specifying one fundamental cosmological pa- 
rameter in a direct way. Nevertheless, we obtained con- 
tours which are significantly tightened for all the param- 
eters, breaking all the degeneracies. This allows us to in- 
vestigate constraints on cosmological parameters without 
introducing any prior from measurements based on stan- 
dard candles hypothesis. In this paper we use Hq — 100ft, 
km/s/Mpc. 

2. Constraints from CMB alone 

2.1. CMB analysis 

During the last decade (since the COBE/DMR results) 
different measurements of CMB anisotropics at differ- 
ent scales have been used to constrain the cosmology of 
our Universe (Hancock et al., 1997 , Lineweaver et al., 
1997, Dodelson and Knox, 2000, Tegmark and Zaldarriaga 
2000ab, Ledour et al., 2000, Jaffe et al., 2001 and others). 
The last experiments (second generation) are sensitive to 
different scales. This gives homogeneous sets of data from 
a few degrees to a few arcmins, which cover two or three 
peaks in the power spectrum with good accuracy. Wc then 
expect the constraints to be stronger, especially if several 
experiments are consistent. 



Our method to derive parameters from CMB data has 
been presented in Bartlett et al., 2000 (BDBL), Douspis 
et al., 2001 (DBBL) and used in Le Dour et al., 2000 
(LDBB)[]. We use the maximisation technique rather than 
the marginalisation obtained by integrating over some pa- 
rameters. Both techniques are equivalent when the prob- 
abilities are Gaussian and the model is linear in the pa- 
rameters. This is far from being the case however when 
dealing with likelihood on the Ce ■ The degeneracies imply 
that likelihood surfaces are sheets in a multidimensional 
space. In the presence of such a complex structure in the 
likelihood space, the marginalisation has the inconvenient 
that the model corresponding to the preferred parame- 
ters may actually lie outside of the actual allowed region! 
We present our results by means of contours in two di- 
mensional parameter space. This allows us to some extent 
to identify complex structure in the multidimensional pa- 
rameter space, that might be hidden in the likelihood on 
one parameter. Goodness of fit has to be evaluated before 
addressing parameter estimation. The technique to derive 
the goodness of fit in our method is detailed in Douspis 
et al. 2001 (DBB). The basics are to use flat-band es- 
timates as published (see table 2 of LDBB) and apply 
our likelihood approximation instead of a non appropri- 
ate x 2 minimisation. Our method remains as simple as a 
X 2 fit but is less biased as shown in DBBL. In this paper, 
we use the most recent results of MAXIMA (Lee et al., 
2001), BOOMERanG (Netterfield et al., 2001) and DASI 
(Halverson et al., 2001). We also use the COBE/DMR 
data (Tegmark and Hamilton, 1997) and refer to the data 
set as CBDM. We also upgraded our grid of models, al- 
lowing closed models and refining the step and range of 
the baryon content (f^ft 2 ), the total density (fltot) and 
the spectral index (n). The parameter space explored is 
given in table 1. 
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Table 1. Parameter space explored: 

Qtot = 1 — f2 re , where £l K is the curvature parameter 

A = A/3 

rjio = (baryon number density) / (photon number density) 

xlO 10 (Note: tl h h 2 = 0.00366r/i ) 
n = pr imeval spe ctral index 
Q = V(5/4^)C 2 



We do not consider reionisation or gravitational waves 
nor neutrino contributions to the matter density f2 m . 
Previous analyses showed that these parameters do not 
play an important role. The present analysis has been 
done with all the information available in the literature. 



1 see http : / / webast . ast . obs-mip ,fr /cosmo/CMB for comple- 
mentary informations 
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Unfortunately, some information is still missing (like the 
window functions); some effects have thus not been taken 
into account but our previous investigation of these ques- 
tions indicates that this lack of information is not critical. 

The likelihood of each model is calculated as described 
in LDBB, and the best model is found by maximising the 
likelihood function over the explored space. The contours 
are defined in the full six-dimensional space with values 
of Alog(£) = 1,4,9 (dashed, in red), corresponding to 
the 1, 2 and 3 a contours for Gaussian distributions when 
projected onto one of the axes. Identically the 1, 2 and 
3 a contours in the two parameters space are defined by 
Alog(£) = 2.3,6.17,11.8 (filled, in blue), for Gaussian 
distribution. Since the likelihood is not actually Gaussian, 
the confidence percentages associated with our contours 
are actually not known; the technique is however standard 
practice. 

We also would like to emphasise the effect of pri- 
ors in the analysis. As it has been noticed in the recent 
literature (Jaffe et al., 2001, Balbi et al., 2000, Lange 
et al., 2000), the final results on the preferred parame- 
ters are very sensitive to priors. Most of the recent work 
set Hq as a "ghost" parameter, using physical densities 
(£l{,/i 2 , ilh 2 , fldmh 2 , ■■■) and deriving Hq by: 
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Then the analysis proceed by putting a prior on H like 
Hq > 45. The representation in one-dimensional plots does 
not allow to know exactly what is the effect of such tech- 



nique. We will show in section 2.2 that such a prior puts 



more severe constraints on fltotai as it is expected. 

The study of the goodness of fit of "best" models is 
detailed in DBB. We apply the same technique on the 
new data and the best model plotted in figure § The 
goodness of fit expressed in terms of generalised \ 2 is 
Xgen = 45.0/43 where 43 = 49 — 6 is the number of degrees 
of freedom corresponding to the number of experimental 
points minus the number of investigated parameters. This 
good value of the goodness of fit allows us to consider 
confidence intervals on cosmological parameters. 



2.2. Constraining the cosmological parameters 

The first Doppler peak is one of the main feature of the 
CMB angular power spectrum. It is clear now that dif- 
ferent experiments indicate both a rise in the power at 
intermediate scale and a fall-off at smaller scales (higher 
£) . The position of the first acoustic peak is strongly re- 
lated to the curvature of the Universe. A statistical anal- 
ysis of the CMB data should then give some strong con- 
straint on the total density (or curvature) of the Universe. 
The confidence intervals on the Sltotai parameter in dif- 
ferent combinations {(Cttotah A), (tltotah H ), (Cl to tal, Vw)) 
are shown as 2-D contour plots in figure The Sltotai — Hq 
plane reveals the stronger constraint that CMB fluctua- 
tions lead to: the contours are almost reduced to a line in 



this plane. The main implication is the exclusion of almost 
all open cosmologies at high confidence level. Preferred 
models are closed (fltotai ~ 1-3) and relative maxima are 
present at lower value (Qtotal ~ 1-05). The models corre- 
sponding to the latter are in good agreement with a flat 
Universe and the inflationary scenario prediction. Within 
the 95% Gaussian confidence level (95% GCL hereafter), 
and allowing Ho to be in then range [20,100] we find: 
0.95 < Sltotai < 1-4. Given our 2-D representation, it is 
easy to see that the closest models have the lowest value 
of Hq (cf. figure ^[). This is a supplementary information 
that the 1-D reduction is unable to show. 

The fourth plot of figures || presents the results in the 
7710 — Sltotai plane. The constraints on 7710 are becoming 
tight with the present day data set compared to analysis of 
first generation experiments (LDBB). They point toward a 
medium value of the baryon content which is well in agree- 
ment with the BBN and light elements determination of 
7710, contrarily to analysis based on first BOOMERanG 
and MAXIMA results (Tegmark & Zaldarriaga, 2000, 
Douspis, 2000, Jaffe et al., 2001). The change in the pre- 
ferred value of flbh 2 is due to the improved measure- 
ments of the second and third peaks. In spring 2000, 
BOOMERanG and MAXIMA experiments measured the 
power at scales corresponding to the second acoustic peak. 
The latter appears to be lower (relative to the first one) 
than expected in a model with a baryon content compati- 
ble with BBN predictions {Q b h 2 ~ 0.02, Buries et al 2001). 
The relative heigh between the two peaks favours high val- 
ues of 7710 corresponding to Vl^ti 2 ~ 0.03. In spring 2001, 
BOOMERanG, DASI and MAXIMA teams are publish- 
ing new analysis showing the possible existence of a third 
acoustic peak with the same height as the second. This dis- 
favours high baryon models which predict a higher third 
peak (see figure 1 of LDBB). The actual value in accor- 
dance with our analysis is 4.0 < 7710 < 7.0 at 95% GCL 
and almost independently of the value of Hq . 

The spectral index of primordial fluctuations, n is also 
a well constrained parameter. We found 0.81 < n < 1.02 
at 95% GCL. This is in agreement with the inflationary 
predictions. 

The remaining parameters are almost not constrained. 
We can see in figure || the different degeneracies between 
parameters. There is no constraint on Ho from the CMB. 
And we can also see that Oa can take almost any value: 
0.0 < tt A < 0.9 at 95% GCL. Another contour plot of 
figure H shows a tight degeneracy between fltotai and Hq . 
There again, it is obvious that the lower Hq, the higher 
fltotai- Imposing Hq > 45 km/s/Mpc results in a tight con- 
straint on fl tot ai ■ 0.9 < fltotai < 1.2 (99% GCL). Finally, 
the second plot of figure || points out the anticorrelation 
of the pair (n, 7710). Increasing n (lowering Q) will increase 
the height of the second peak relative to the first. At the 
opposite, increasing 7710 will decrease the second peak. 

It is clear that the CMB alone does not directly con- 
strain very tightly Sl to tai without further restriction (or 
prior) on Hq. Nevertheless it is important to keep in 
mind that the CMB impose a very narrow region in the 
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Fig. 2. Different contour plots from the analysis of set CBDM. The dashed red line define the 68, 95 and 99% GCL 
when projected on the axis. The blue filled contours are the corresponding confidence intervals in 2 dimension. 



^totai ~ H a plane. Other parameters are not so well con- 
strained. However, the baryon content is now becoming 
rather well constrained, in agreement with nucleosynthesis 
and light elements abundances and the primordial index 
is slightly below 1, close to the value predicted by infla- 
tion. It is therefore a fundamental success of Cosmology 
that the CMB alone does allow such tight constraints. 
However, it is not possible yet to infer robust constraint 
on other fundamental parameters like Hq, Q\ or f2 m . It 
is therefore vital to have other means to measure these 
fundamental parameters. 

3. Combined analysis 

Given the fact that CMB implies degeneracies between 
the fundamental cosmological parameters Hq, fl\ and Q. m , 
it is interesting to look for combination with constraint 
that do not measure directly any of these parameters. 
Primordial nucleosynthesis is one of this possibility; we 
have checked that this does not change much the con- 
straints established in the above section. An other possi- 
bility that we investigate below is to use the baryon frac- 
tion inferred from cluster observations. 

3.1. Baryon fraction in clusters 

Clusters provide fundamental observations for cosmology. 
Historically, they provided the first evidence for the ex- 
istence of dark matter. Furthermore they now provide us 
with the best evidence that this dark matter is non bary- 
onic. Indeed, clusters are the only structure for which the 
total mass, baryon mass and stellar mass can be evalu- 
ated simultaneously in a reliable way. Cosmological appli- 



cations of clusters have been revived with new fundamen- 
tal tests which provide global constraint on the cosmo- 
logical parameters: the baryon fraction in clusters trace 
reliably the ratio between baryons and the total mass in 
the Universe: 



where T has been estimated from numerical simulations 
and is close to 0.92 in the outer part of clusters (Frenk et 
al., 1999). Provided that fif, can be estimated from primor- 
dial nucleosynthesis and the observed abundance of light 
elements, this argument can be used to estimate a robust 
upper limit on Slo (White et al., 1993). In this paper we 
combine CMB and the sole baryon fraction constraints. 
This is interesting because the baryon fraction does not 
constrain by itself any of the cosmological parameters, so 
that in principle it could leave unchanged the conclusions 
obtained from the CMB alone. 

3.2. Combining constraints 

There are some dispersions among values published for 
the gas fraction in clusters, which mainly reflect the dif- 
ference in the mass estimator used. Using mass estima- 
tors based on numerical simulations a median baryon frac- 
tion of 0.048ft- 3 / 2 + 0.014 (±10%) was recently found by 
Roussel et al., 2000, consistently with previous investiga- 
tions (Arnaud & Evrard, 1999). More recently the stan- 
dard value has been revisited by Sadat and Blanchard, 
2001, showing that the actual baryon fraction could actu- 
ally be lower. They concluded that a baryon fraction as 
low as 0.031ft.- 3 / 2 + 0.012 (±10%) is still consistent with 
the X-ray data on clusters. 
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Fig. 3. Combined analysis of CMB and high baryon frac- 
tion. See figure for the definition of contours. 



Fig. 4. Combined analysis of CMB and low baryon frac- 
tion. See figure ||for the definition of contours. 



3.2.1. The high baryon fraction case 

As we have mentioned previously, we have chosen to in- 
vestigate the combination of the baryon fraction /& with 
the CMB data. The baryon fraction is an interesting case 
because its value can be measured in a rather direct way 
through observations of X-ray clusters provided the as- 
sumption that all the baryons are actually seen (assum- 
ing no large quantity of dark baryons) . It is worth notic- 
ing that without an additional constraint on fi;,, there is 
no trivial reason why the additional information on 
should further constrain the various cosmological param- 
eters. However, it happens that the additional constraint 
from the baryon fraction restricts significantly the various 
contours (see figure ||). This is very clear from the H — i]i 
plane: the area of the contours is much reduced, leading 
to a preferred model which is around Hq ~ 60 km/s/Mpc 
and 7710 ~ 5. The other two-dimensional plots also reveal 
that all the contours are significantly reduced: actually 
most of the one sigma contours are of the order of the 
grid size so that best model and one sigma ranges should 
be used with caution. It is highly interesting that the cos- 
mological model is now very well constrained: il\ ~ 0.6, 
fi tot - 1. (implying VL m ~ 0.4, H ~ 60 km/s/Mpc, 
7710 ~ 5, n ~ 0.85). As we mentioned the allowed range 
have to be interpreted with caution, but we notice that the 
preferred value of the Hubble constant is rather low and 
that higher values like 80 km/s/Mpc lie uncomfortably 
outside of the preferred region. 

3.2.2. The low baryon fraction case 

As mentioned above, a baryon fraction significantly lower 
than current estimates is still consistent with, and actually 



preferred by existing data on X-ray clusters. Using this low 
value for the baryon fraction leads similarly to much re- 
stricted contours. However the latter differ from the previ- 
ous ones (see figure ||) : the preferred value for 7710 is again 
close to 5., the preferred value of the Hubble constant is 
now very low (Hq ~ 40 km/s/Mpc). Although such low 
values have been argued from time to time (Bartlett et 
al., 1995), it is worth noticing that values as Hq ~ 60 
km/s/Mpc are still at the edge of the 3 sigma contour and 
for this reason cannot be entirely ruled out. The preferred 
model is tt\ ~ 0.3, Cttot ~ 1-1 (implying and VL m ~ 0.8, 
H ~ 40 km/s/Mpc, r? 10 ~ 5., n ~ 0.85). 

Finally, we also investigated the additional constraint 
that can be inferred from the evolution of the abundance 
of X-ray clusters: the evolution of the abundance of X-ray 
clusters is known to be a powerful constraint on the matter 
density of the universe (Oukbir and Blanchard, 1992). The 
measurement of the temperature of a flux-limited sam- 
ple of X-ray clusters has allowed the determination of the 
abundance of clusters at z ~ 0.3 — 0.4 (Henry, 1997; Henry, 
2000) . Using the updated estimates of local abundance of 
X-ray clusters a significant level of evolution is found, in- 
dicating a high density Universe (Blanchard et al., 2000). 

The constraint provided by the addition of the evolu- 
tion of cluster abundance is essentially a constraint on 
Q m - Therefore it is rather natural that this additional 
constraint by itself will not change much the previous 
contours as the value found by Blanchard et al., 2000, 
Q m ~ 0.8, is almost identical to what is provided by the 
previous analysis. We can anticipate that using analysis 
leading to lower O m with the higher baryon fraction will 
not change much the contours presented in the previous 
section. 
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4. Conclusions 

The new CMB data provide impressive evidence in favour 
of the standard inflationary picture for structure forma- 
tion with non baryonic dark matter. They also brilliantly 
confirm previous evidence for a nearly flat Universe. In this 
paper, we have actually shown that the data provide a very 
narrow allowed region in the Hq — f2 tot plane. However, it 
is somewhat frustrating that despite the high quality of 
the data, degeneracies among the fundamental cosmolog- 
ical parameters (D, t oti^\ and Hq) cannot be broken and 
that constraints on other parameters are relatively weak. 
It should also be realised that given the dependency on 
the power spectrum shape and the possible contribution 
at some level from topological defects, it will be extremely 
difficult to obtain entirely reliable constraints on these pa- 
rameters from CMB alone. Actually, even with the preci- 
sion anticipated for the Planck mission, the degeneracy 
problem may not be solved. It is therefore vital to develop 
new tests that provide additional constraints. Supernova 
Hubble diagram provide such a test which has be widely 
advertised. However, the previous use of the Hubble di- 
agram has always been polluted by systematic effects, 
and therefore we believe that a coherent picture will be 
firmly established only from further evidence. Clusters of- 
fer such additional possibility. We have therefore investi- 
gated the complementary information which come from 
clusters. The most interesting constraint has been found, 
rather surprisingly, to be the baryon fraction. Indeed with 
this single additional constraint we found that nearly all 
the cosmological parameters are well specified. Using a 
classical value for the baryon fraction of ~ 15% for h = 0.5 
we obtain as the best model: Q\ ~ 0.6 , Q to t ~ 1 ■, Hq ~ 60 
km/s/Mpc, 7710 ~ 5., n ~ 0.85. Using a more recent re- 
sult from Sadat and Blanchard, 2001, who infer a lower 
baryon fraction when several systematics are corrected for, 
we found that the best model is fl\ ~ 0.3, Cltat ~ 1.1 > 
Ho ~ 40 km/s/Mpc, 7710 ~ 5., n ~ 0.85 implying a high 
density parameter Q m ~ 0.8 consistent with the deter- 
mination from cluster abundance evolution by Blanchard 
et al., 2000. Although the Einstein-de Sitter model is at 
the edge of the 3 sigma contour, it is probably premature 
to rule it out on this basis. It is remarkable that in both 
analyses we found that the preferred primordial baryon 
content is narrowly constrained and consistent with pri- 
mordial nucleosynthesis. Identically the primordial index 
is found to be of the order of 0.85 in both cases. 
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